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ABSTRACT 

We present time-resolved spectroscopy and photometry of the double-lined eclipsing 
cataclysmic variables AC Cnc and V363 Aur (= Lanning 10). There is evidence of 
irradiation on the inner hemisphere of the secondary star in both systems, which we 
correct for using a model that reproduces the observations remarkably well. We find 
the radial velocity of the secondary star in AC Cnc to be Kji = 176 ± 3 kms~^ 
and its rotational velocity to be wsini — 135 ± 3 kms~^. From these parameters 
we obtain masses of Mi = 0.76 ± 0.03 Mq for the white dwarf primary and M2 = 
0.77 ± 0.05 Mq for the K2 ± IV secondary star, giving a mass ratio of g = 1.02 ± 
0.04. We measure the radial and rotational velocites of the G7 ± 2V secondary star in 
V363 Aur to be Kn = 168 ± 5 kms~^ and wsin j = 143 ± 5 kms""'^ respectively. The 
component masses of V363 Aur are Mi = 0.90 ± 0.06 Mq and M2 = 1.06 ± 0.11 Mq, 
giving a mass ratio of g = 1.17 ± 0.07. The mass ratios for AC Cnc and V363 Aur fall 
within the theoretical limits for dynamically and thermally stable mass transfer. Both 
systems are similar to the SW Sex stars, exhibiting single-peaked emission lines with 
transient absorption features, high-velocity S-wave components and phase-offsets in 
their radial velocity curves. The Balmer lines in V363 Aur show a rapid increase in 
flux around phase followed by a rapid decrease, which we attribute to the eclipse 
of an optically thick region at the centre of the disc. This model could also account 
for the behaviour of other SW Sex stars where the Balmer lines show only a shallow 
eclipse compared to the continuum. 

Key words: accretion, accretion discs - binaries: eclipsing - binaries: spectroscopic - 
stars: individual: AC Cnc - stars: individual: V363 Aur - novae, cataclysmic variables. 



1 INTRODUCTION 

Cataclysmic variables (CVs) are close binary stars consist- 
ing of a red dwarf secondary transferring material onto a 
white dwarf primary via an accretion disc or magnetic ac- 
cretion stream. AC Cnc and V363 Aur are both examples 
of nova-likes (NLs), defined as CVs which have never been 
observed to undergo nova or dwarf-nova type outbursts. See 
[Warner (1995a) [for a comprehensive review of CVs. 
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A knowledge of the masses of the component stars in 
CVs is fundamental to our understanding of the origin, evo- 
lution and behaviour of these systems. Population synthesis 
models and the disrupted magnetic braking model of CV 
evolution can be observationally tested only if the number 
of reliably known CV masses increases. One of the most 
reliable ways to measure the masses of CVs is to use the 
radial velocity and the rotational broadening of the sec- 
ondary star in eclipsing systems; the radial velocity of the 
disc emission lines may not represent the white dwarf's or- 
bital motion. At present, reliable masses are known for only 
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~20 CVs, partly due to the difficulties in measurement (see 
ISmith fc Dhillon 19981 for a review). 

AC Cnc was classified as an eclipsing variable 
star by [Kurochkin fc Shugarov (1980)| with an or- 
bital period of 7^2 hours. [Shugarov (1981)| suggested 
that AC Cnc is a NL based on UBV colours, and 
the CV nature of AC Cnc was confirmed spectro- 
scopically by [Okazaki, Kitamura fc Yamasaki (1982)| 
through broad H and He emission lines and the 
later study of [Yamasaki, Okazaki fc Kitamura (1983)| 
[Downes (1982)[ discovered secondary star features in 
the spectra leading to the first mass determination by 
[Schlegel, Kaitchuck fc Honeycutt (1981)1 who found Afi = 
0.82 ± 0.13 Mq and M2 = 1.02 ± 0.14 Mq from the radial 
velocities of the primary and secondary components. 

V363 Aur {— Lanning 10) was discovered by 
[Lanning (1973)[ as a UV-bright source and later 
found to be a CV through broad Balmer and Hell 
emission by [Margon fc Downes (1981)[ and its typ - 
ical C V energy distribution I jSzkody fc Crosa 1981| l . 
[Home, Lanning fc Comer (1982)1 obtained the first spectro- 
scopic and photometric data, finding that V363 Aur 
is an eclipsing system with a period of 7.7 hours. 
[Schlegel, Honeycutt fc Kaitchuck (1986)] calculated the 
component masses of V363 Aur to be Mi — 0.86 ± 0.08 
Mq and M2 = 0.77 ± 0.04 Mq from the radial velocities of 
the HeH A4686A emission line and the G-band absorption. 

The existing component masses of AC Cnc and 
V363 Aur use emission line radial velocity curves that ex- 
hibit phase shifts, suggesting that they could be unreliable. 
In addition to this, the mass ratio of AC Cnc found by 
[Schlegel, Kaitchuck fc Honeycutt (i984j| of q = 1.24 ± 0.08 
is the highest known of any CV and is very close to the 
upper limit of mass transfer stability computed by the mod- 
els of [Politano (1996)] In this study, we present photometry 
and spectroscopy of AC Cnc and V363 Aur to calculate new 
masses using the secondary star properties alone. 



The spectra and images were reduced using stan- 
dard procedures (e.g. [Dhillon, Jones fc Marsh 1994| 
[Thoroughgood et al. 200T] . Comparison arc spectra were 
taken every 40-50 min to calibrate instrumental flexure. 
The arcs were fitted with a sixth-order polynomial in 
blue and a fourth-order polynomial in red with rms 
scatters of better than O.OlA. The photometry data were 
corrected for the effects of atmospheric extinction by 
subtracting the magnitude of a nearby comparison star 
(AC Cnc-9 and V363 Aur-3; [Henden fc Honeycutt 1995t 
and using values obtained by the CAMC telescope 
C Helmer fc Morrison 1985ll . The absolute photometry is ac- 
curate to approximately ±0.5 mjy; the relative photometry 
±0.03 mag. Slit losses were then corrected for by dividing 
each AC Cnc and V363 Aur spectrum by the ratio of the 
flux in the spectrum (summed over the whole spectral 
range) to the corresponding photometric flux. 



3 RESULTS 

3.1 Ephemeris 

We derived new ephemerides for AC Cnc and V363 Aur, 
which are used to phase all data presented in this paper. The 
times of mid-eclipse were determined by fitting a parabola 
to the eclipse minima in the JKT data. 

In the case of AC Cnc, a least-squares fit to the 29 
eclipse timings listed in Table [5] (a) yields the ephemeris: 



rmid-cciipso = HJD 2 451924.5368 ± 0.30047747 E 
± 0.0006 ± 0.00000004. 



(1) 



A least-squares fit to the 17 eclipse timings of V363 Aur 
listed in Table [21(b) gives the ephemeris: 



rmid- 



cclipsc 



HJD 2 444 557.9514 ± 0.32124187 E 
± 0.0016 ± 0.00000008. 



(2) 



We see no evidence for any systematic variation in the 0-C 
values shown in Table [2| in either AC Cnc or V363 Aur. 



2 OBSERVATIONS AND REDUCTION 

On the nights of 2001 January 9-14 we obtained blue and 
red spectra of AC Cnc and V363 Aur with the 2.5-m Isaac 
Newton Telescope (INT) ± IDS spectrometer on La Palma. 
The blue setup comprised of the 235-mm camera with the 
R1200B grating and the EEVIO CCD chip, which gave a 
wavelength coverage of approximately 4490-5580A at 0.95- 
A (57 kms~^) resolution. In the red we used the 500-mm 
camera with the R1200Y grating and the TEK5 CCD chip 
resulting in a wavelength range of 6320-6720A at a resolu- 
tion of 0.8-A (36 kms~^). Simultaneous photometry in the 
Johnson-Cousins B and R bands was recorded with the 1- 
m Jacobus Kapteyn Telescope (JKT) using the SITe2 CCD 
chip. Full phase coverage was achieved for both objects - a 
full journal of observations is given in Table^ including the 
exposure times used. 

We also collected 19 spectral type templates ranging 
from G5V-M2V, telluric stars to remove atmospheric fea- 
tures and flux standards on both the INT and JKT. Seeing 
varied between 1.0 and 1.5 arcsec over the observing run and 
conditions were photometric on all nights except for January 
10 when some patchy cloud was present. 



3.2 Average Spectrum 

The average spectra of AC Cnc and V363 Aur are shown in 
Fig-H fl'iid in Table|2|we list fluxes, equivalent widths (EW) 
and velocity widths of the most prominent lines measured 
from the average spectra. 

Both systems show broad, symmetric, single-peaked 
Balmer emission lines instead of the double-peaked pro- 
files one would expect from a high inclination accretion disc, 
much like other nova-like systems (e.g. lWarner 1995at . The 
Hel lines, however, are broad and double-peaked in na- 
ture. The line strength of Hell A4686A is much stronger 
in V363 Aur than AC Cnc and even more dominant than 
H/3 emission. Another high excitation feature, the CIII/NIII 
AA4640-4650A blend, is only present in V363 Aur and is very 
broad. The emission lines are characteristic of the SW Sex 
stars (e.g. [Dhillon, Marsh fc Jones 1997] , but unlike others 
in this sub-class, these systems show clear secondary star 
features (no doubt due to their longer periods and hence 
earlier type secondaries). Both AC Cnc and V363 Aur show 
the absorption features of the neutral metals Cal, Fel and 
Mgl, even in the average spectrum shown in Fig. which 
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Table 1. Journal of observations for AC Cnc and V363 Aur. The epochs are calculated using the new ephemerides presented in this 
paper (equations 111 and l2t . 



UT Date 



Object 



INT 
setup 



No. of 
spectra 



Exposure 
time (s) 



JKT 

filter 



No. of 
images 



Exposure 
time (s) 



UT 
start 



UT 
end 



Epoch 
start 



Epoch 
end 



2001 Jan 09 


V363 Aur 


Red 


61 


2001 Jan 10 


V363 Aur 


Blue 


39 


2001 Jan 11 


V363 Aur 


Blue 


70 


2001 Jan 12 


V363 Aur 


Red 


64 


2001 Jan 12 


AC Cnc 


Red 


33 


2001 Jan 13 


AC Cnc 


Red 


96 


2001 Jan 14 


AC Cnc 


Blue 


90 



300 


R 


323 


400 


B 


272 


400 


B 


474 


300 


R 


348 


300 


R 


179 


300 


R 


530 


300 


B 


521 



30 


22:35 


03:55 


22915.74 


22916.43 


30 


20:07 


04:36 


22918.53 


22919.63 


30 


19:33 


04:06 


22921.57 


22922.67 


30 


22:32 


04:11 


22925.97 


22926.79 


30 


04:26 


07:21 


-6.14 


-5.74 


30 


22:39 


07:04 


-3.31 


-2.45 


30 


22:27 


06:33 


-0.32 


0.80 



has not been corrected for orbital motion. The secondary 
star features appear to be stronger relative to the contin- 
uum in AC Cnc than in V363 Aur. The weak feature at 
A6614A in the V363 Aur red spectrum is an interstellar 
absorption line. 



3.3 Light Curves 

The top panels of Fig.|2ta) and I3b) show the B and i?-band 
JKT light curves. The remaining panels show emission-line 
light curves, which were produced by subtracting a polyno- 
mial fit to the continuum and summing the residual flux. All 
light curves are plotted as a function of phase following our 
new ephemerides. 

The B and i?-band JKT light curves of AC Cnc show 
deep, symmetrical primary eclipses with out-of-eclipse mag- 
nitudes of 14.30 ± 0.05 mag in B and 14.00 ± 0.05 in R. 
The primary eclipse depth is 1.8 mag in B and 0.73 mag 
in R. We measured the phase half-width of eclipse at the 
out-of-eclipse level (A^) by timing the first and last con- 
tacts of the B and i?-band photometry eclipses and divid- 
ing by two. Our average value of A^ = 0.09 ± 0.01 phases 
is smaller than, but consistent with, the value of 0.108 ± 
0.008 quoted in |Harrop-Allin fc Warner (1996)| Ellipsoidal 
modulation of the red star is clearly present, although fiaring 
around phases 0.3-0.4 contaminates the effect in B. There 
is also evidence for a secondary eclipse at phase 0.5. We see 
no evidence for a bright-spot in the light curves but fiicker- 
ing is present, particularly just after primary eclipse. A no- 
table feature of the B and i?-band JKT Ught curves is the 
U-shaped eclipse minima, in contrast to the V-shaped min- 
ima seen in many SW Sex systems ( |Knigge et al. 2000| l . The 
eclipses of the Balmer lines show the usual V-shape, but the 
Hell line has a U-shaped eclipse minimum and is completely 
eclipsed, suggesting an origin close to the white dwarf. The 
Ha fiux increases markedly after eclipse before slowly de- 
clining - there is also the suggestion of a sharp decrease 
in flux around phase 0.5. This secondary eclipse is possibly 
also seen in the Hel line, although the primary eclipse here is 
much broader and shallower. The H/3 flux seems more erratic 
in behaviour, closely resembling the higher-excitation HeH 
emission line. Note that when the Balmer flux increases, the 
fiickering in the JKT Ught curves is more prominent. 

The V363 Aur JKT light curves in both the B and 
7?-bands are deep and symmetrical with V-shaped eclipse 
minima, much like the SW Sex systems ( |Knigge et al. 2000^ . 
The B-band out-of-eclipse magnitude is 14.50 ± 0.10, 
and the eclipse depth is 0.88 mag; in the 7?-band, the 
out-of-eclipse magnitude is 13.65 ± 0.05 with an eclipse 



depth of 0.63 mag. Our measured phase half-width of 
eclipse at the out-of-eclipse level, Ai^ = 0.078 ± 0.005 
(the average from the B and R-h&nd photometry eclipses) 
is lower than the value of 0.120 ± 0.010 quoted by 
[Harrop-AUin fc Warner (1996)| There is evidence of either 
a shallow secondary eclipse or orbital modulation in the 
7?-band Ught curve but not so in B. One of the most no- 
table features of the UgUt curves is the high level of fiick- 
ering present. The red emission-line Ught curves of Ha and 
Hel are similar in that they both show a maximum fiux 
around phase 0.7. Perhaps the most interesting feature of 
the emission-line lightcurves is in the primary eclipse of the 
Balmer lines; the fiux seems to drop entering eclipse but at 
phase there in a sharp increase in Balmer line emission 
followed by a rapid decrease. This is particularly prominent 
in H/9, but also seems to be present in the Ha line and pos- 
sibly the Hel line. The effect is definitely not present in the 
high excitation HeH and CHI/NHI complex. 



3.4 Trailed spectrum &: Doppler Tomography 

We subtracted polynomial fits to the continuum from the 
spectra and then rebinned the spectra onto a constant 
velocity-interval scale centred on the rest wavelength of the 
Ha, H/3, Hel A6678A and Hell A4686A lines. The data were 
then phase-binned into 50 bins, all of which were filled ex- 
cept for 1 empty bin at phase 0.31-0.33 in the blue V363 Aur 
data set. The trailed spectra of the lines are shown in the 
upper panels of Fig. |3] (a) and (b). We then constructed 
Doppler tomograms from the trailed spectra, a technique 
which maps the velocity-space distribution of the emission 
lines (e.g. IMarsh 2001t . The Doppler maps are shown in 
Fig. 2] and trailed spectra reconstructed from the maps are 
presented in the lower panels of Fig. |3 (a) and (b) . 

The trailed spectra of AC Cnc show two clearly-defined 
components. The flrst is a high-amplitude S-wave with a 
semi-amplitude of ~500kms~^, which crosses zero veloc- 
ity from red to blue around phase 0.15. The emission is 
particularly noticable in Ha but can be seen in the other 
lines. This component appears in the Doppler map super- 
imposed upon a ring of emission characteristic of an ac- 
cretion disc at Rd ~ 0.4 — 0.5Li. This is what one would 
expect from the bright-spot, although it appears sligUtly 
downstream from where the computed gas stream joins the 
accretion disc, a trait also seen in other CVs (e.g. WZ Sge; 
[Spruit fc Rutten 1998| l. The second, lower-velocity compo- 
nent seen in the Balmer lines, seemingly in anti-phase with 
the higher-velocity component, shows up clearly in the 
Doppler maps as originating on the inner hemisphere of the 
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Figure 1. Average spectra of AC Cnc (left) and V363 Aur (right). 



w 



Table 3. Fluxes and widths of prominent lines in AC Cnc and V363 Aur, measured from the average spectra. In the case of V363 Aur, 
Hell A4686A and CIII/NIII AA4640-4650A are blended, so separate values of the flux and EW arc given (determined from a double- 
Gaussian fit) as well as the combined flux of the two. 

(a) AC Cnc 



Line 



Flux 

X 10"" 

(ergscm~^ s^^) 



EW 

(A) 



FWHM 
(kms-i) 



FWZI 
(kms-i) 





7.57 ±0.01 


16.55 ±0.03 


1100 ±100 


2900 ± 500 




4.63 ±0.01 


8.57 ±0.03 


1250 ± 100 


2800 ± 300 


A6678A 


0.94 ± 0.01 


2.17 ±0.02 


1100 ±100 


1900 ± 500 


I A4686A 


1.36 ±0.01 


2.32 ± 0.02 


1600 ± 100 


2300 ± 300 



(b) V363 Aur 



Ha 8.07 ±0.02 11.82 ±0.02 1150 ±100 3100 ± 500 

H/3 3.60 ± 0.02 4.88 ± 0.02 1250 ± 100 2900 ± 300 

HeIA6678A 0.86 ± 0.01 1.31 ± 0.02 1100 ±100 2300 ± 500 

HeIIA4686A 4.84 ± 0.02 6.19 ±0.03 1150 ±100 3400 ± 500 

CIII/Nin AA4640-4650A 2.05 ± 0.03 2.58 ± 0.03 1950 ± 100 4900 ± 500 

Hell + Cni/NIII 7.35 ± 0.02 9.37 ± 0.03 



secondary star. There is possibly another component visi- 
ble in the Hel A6678A line, which could be interpreted as 
the signature of a faint double-peaked accretion disc. All of 
the low-excitation lines exhibit the rotational disturbance 
expected from a high inclination accretion disc, but inter- 
estingly this is not seen in the Hell A4686A line. These 
emission features are very similar to those seen in the nova- 
like SW Sex dPhillon, Marsh fc Jones 1997| l. 

The high-amplitude S-wave and the low-amplitude 
emission are also seen in V363 Aur, although the anti-phase 
sinusoid of the secondary is less well pronounced. The high- 
amplitude Balmer emission in the Doppler maps appears as 
a bright -spot downstream from where the gas stream meets 
the disc, as well as from the opposite side of the disc. There 
is clearly Balmer emission from the inner hemisphere of the 
secondary which is likely to be caused by irradiation by the 
accretion regions or white dwarf. The Hel trailed spectrum 



gives the visual impression of an absorption feature (moving 
at ~ Kw) super-imposed on a high-velocity S-wave. The 
strange increase in Balmer line and Hel A6678A emission 
at phase zero seen in the lightcurves is clearly visible in the 
trailed spectra. Possible explanations for this are given in 
Section 0] The high-excitation HeH A4686A line is single- 
peaked throughout the orbit and the emission lies on the 
white dwarf centre of mass in the Doppler map. 



3.5 Radial velocity of the white dwarf 

We measured the radial velocities of the emission lines in 
AC Cnc and V363 Aur by applying the double-Gaussian 
method of [Schneider fc Young (1980)| since this technique is 
sensitive mainly to the line wings and should therefore reflect 
the motion of the white dwarf with the highest reliability. 
We used Gaussians of widths 200, 300 and 400 kms~^ and 
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Figure 2. Lightcurves of (a) AC Cnc and (b) V363 Aur. Top panels: R and B-band JKT lightcurves. Middle and bottom panels; 
Emission-line lightcurves of Ha, H/3, CIII/NIII AA4640-4650A + Hell A4686A and Hel A6678A. The i?-band V363 Aur data are from 
2 nights; open circles represent data from 12/01/01 and closed circles represent data from 09/01/01. Note that the two nights of data 
are split by orbital phase, as well as by night of observation. 
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Figure 3. (a) Trailed spectra of AC Cnc (upper panel) & data computed from the Doppler maps (lower panel), (b) Trailed spectra 
of V363 Aur (upper panel) & data computed from the Doppler maps (lower panel). A comparison of the upper and lower panels gives 
some indication of the quality of the fit. The gaps in the computed data correspond to eclipse spectra (selected using the light curves of 
Fig.|3J, which have been omitted from the fit as Doppler tomography cannot properly account for these phases. Note that the data have 
been folded to show more than one binary cycle. 
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Figure 4. Doppler maps of emission lines in AC Cnc (left) and V363 Aur (right). The three crosses marked on each Doppler map 
represent the centres of mass of the secondary (upper cross), the system (middle cross) and the primary (lower cross). These crosses, 
the Roche lobe of the secondary star and the predicted trajectory of the gas stream have been plotted using the X^-corrected system 
parameters summarised in Table Isl The series of circles along the gas stream mark the distance from the white dwarf at intervals of 
O.lLi, ranging from l.OLi at the red star to O.SLi. 



varied their separation from 200 to 2500 kms ^. We then 
fitted 



V = ')-Ksm[2%{(i}-(l}o)] 



(3) 



to each set of measurements, where V is the radial velocity, 
K the semi-amplitude, (j) the orbital phase, and 0o is the 
phase at which the radial velocity curve crosses from red to 
blue. 

Examples of the radial velocity curves for AC Cnc 
and V363 Aur are shown in Fig. |S] The most striking 
feature of all of the radial velocity curves are the phase 
shifts, where the spectroscopic conjunction of each line 
occurs after photometric eclipse. This phase shift implies 
an emission line source trailing the accretion disc, such 
as a bright spot, and is a common feature of SW Sex 
stars (e.g. D W UMa, [SB aftcr, Hcssman fc Zha ng 1988] 
V1315 Aql, IDhillon, Marsh fc Jones 1991. SW Sex, 
[Dhillon, Marsh fc Jones 1997^ . There is clear evidence of 
rotational disturbance in the Balmer lines of AC Cnc, 
where the radial velocities measured just prior to eclipse 
are skewed to the red, and those measured after eclipse are 
skewed to the blue. This confirms the detection of a similar 
feature in the trailed spectra, and indicates that at least 
some of the emission must originate in the disc. We tried to 
measure white dwarf radial velocity {Kw) values from the 
emission lines in AC Cnc and V363 Aur using a diagnostic 
diagram (e.g. [Shafter, Szkody fc Thorstensen 1986| l and 
a hght centres diagram (e.g. IMarsh 198811 but with no 
success. This is not surprising given that the Doppler maps 
(Fig. |1J show that the accretion disc does not dominate the 



emission in these systems. We conclude that the emission 
lines of AC Cnc and V363 Aur are unreliable indicators of 
the white dwarf radial velocity due to the phase shifts. 

3.6 Radial velocity of the secondary star 

The secondary star in both AC Cnc and V363 Aur 
is clearly visible in Fig. Q through absorption lines of 
Mgl, Pel and Cal. We compared regions of the spectra 
rich in absorption lines with a number of red dwarfs of 
spectral types G5V-M2V. A technique known as skew 
mapping was used to enhance the secondary features 
and obtain a radial velocity {Kr) measurement (see 
[Vande Putte et al. (2003) | for a detailed cr itique of the 
method and [Smith, Dhillon fc Marsh (1998)| for a success- 
ful application to BT Mon). 

The first step was to shift the spectral type template 
stars to correct for their radial velocities. We then normal- 
ized each spectrum by dividing by a first-order polynomial 
fit, and then subtracting a higher order fit to the contin- 
uum. This ensures that line strength is preserved along the 
spectrum. The AC Cnc and V363 Aur spectra were normal- 
ized in the same way. The template spectra were artificially 
broadened to account for the orbital smearing of the CV 
spectra due to their exposure times (texp) using the formula 

V = ^--v'^'^^i^- (4) 



(e.g. IWatson fc Dhillon 200l( . and then by the rotational 
velocity of the secondary (v sin i) . Estimated values of Kr 
and V sin i were used in the first instance, before iterating 
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Figure 5. Radial velocity curves of the emission lines in AC Cnc using Gaussian widths of 200 kms~^ and a separation of 1200 
knis~^ (left), and of V363 Aur using Gaussian widths of 300 kms~^ and a separation of 1200 kms~^(right). We omitted the points 
around primary eclipse (open circles) during the fitting procedure as these measurements are affected by the rotational disturbance. The 
horizontal dashed lines represent the 7 velocities given by the sinusoidal fits. 



to find the best-fitting values given in Section 13.101 Re- 
gions of the spectrum devoid of emission lines were then 
cross-correlated with each of the templates yielding a time 
series of cross-correlation functions (CCFs) for each tem- 
plate star. To produce the skew maps, these CCFs were 
back-projected in the same way as time-resolved spectra 
in standard Doppler tomography IMarsh fc Home 198811 . If 
there is a detectable secondary star, we expect a peak at 
{OjKh) in the skew map. This can be repeated for each of 
the templates, and the final skew map is the one that gives 
the strongest peak. 

The AC Cnc skew maps show well-defined peaks at 
Ky ~ 186 kms~^ - the skew map for the K2V template 
is shown in Fig. [5] together with the trailed CCFs and the 
regions used for the cross-correlation can be seen in Fig. |7| 
A systemic velocity of 7 = 40 kms~^ was applied in order 
to shift the skew map peaks onto the K^ = axis. The 
Kr value varies little with 7 in practice, as Kx <^ Ky in 
the back-projections (e.g. [Smith, Dhillon fc Marsh 1998| l. 
We adopt 7 = 40 ± 5 kms~ as the systemic velocity of 
the AC Cnc, in contrast to values of 122 kms~^ and 107 
kms~^ found by [Schlegel, Kaitchuck fc Honeycutt (1984)] 
using Balmer emission lines. Our adopted Kr value of 186 
± 3 kms~^ is derived from the best-fitting template (K2V), 
with the error incorporating the spread of values obtained 
by using different templates (see Table |3|l. The uncertainty 
also refiects the scatter in the radial velocity fits shown in 
Fig. |S] Note the remarkable agreement in the {Ky,Kx) val- 
ues obtained from the red and blue data sets in Table|l] The 
small scatter in Ky assures us that Kr is robust, and the 
small scatter in K^ around zero indicates that our assumed 
7 is correct. 

The final V363 Aur skew maps for the G7V template 
(blue) and G5V template (red) and trailed CCFs are shown 
in Fig. |S1 with the regions used for the cross-correlation 
marked in Fig. |7] For V363 Aur, the systemic velocity 
was less simple to determine, as the blue skew map sug- 
gested 7 = -10 kms~^ and the red skew map gave 7 = 



10 kms^^. We can find no explanation for this discrep- 
ancy, so adopt a systemic velocity of 7 = ± 10 kms^^. 
[Schlegel, Honeycutt fc Kaitchuck (1986) ] find systemic ve- 
locities ranging from 10 ± 3 kms~ for the HeH A4686A 
line-wings to 35 ± 3 kms~^ for the G-band absorption line. 
However, the peaks in the skew maps appear consistently 
at Ky PS 184 kms"^. Our adopted value of Kr = 184 ± 
5 kms^^ once again acknowledges the uncertainty in us- 
ing different templates and the scatter in the radial velocity 
curves. 



3.7 Rotational velocity and spectral type of the 
secondary star 

In order to maximise the strength of the secondary fea- 
tures, we averaged the orbitally-corrected eclipse spectra 
of AC Cnc and V363 Aur. The spectral-type templates 
were broadened for smearing due to orbital motion as be- 
fore and rotationally broadened by a range of velocities (50- 
240 kms~^). We then ran an optimal subtraction routine, 
which subtracts a constant times the normalized template 
spectrum from the normalized, orbitally-corrected CV spec- 
trum, adjusting the constant to minimize the residual scat- 
ter between the spectra. The scatter is measured by carrying 
out the subtraction and then computing the x^ between the 
residual spectrum and a smoothed version of itself. By find- 
ing the value of rotational broadening that minimizes the 
X^, we obtain an estimate of both i;sini and the spectral 
type of the secondary star. Note that the vsini values of 
the template stars are much lower than the instrumental 
resolution, so do not affect our measurements of vsini for 
the secondary star. 

The value of v sin i obtained using this method varies 
depending on the spectral type template, the wavelength 
region for optimal subtraction, the amount of smoothing 
of the residual spectrum in the calculation of x^ a-nd the 
value of the limb-darkening coefficient used in the broaden- 
ing procedure. The values of v sin i found from the G and K 
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Figure 6. Left; Skcwmaps and trailed CCFs of AC Cnc, measured using a K2V template. Right: Skewmaps and trailed CCFs of V363 
Aur, measured using a G7V template in the blue and a G5V template in the red. 



templates in the red and blue wavelength ranges, calculated 
using a limb-darkening coefhcient of 0.5 and smoothed us- 
ing a Gaussian of FWHM — 15km s~^, are listed in TableEl 
together with the minimum x^ . The optimal subtraction 
technique also tells us the value of the constant by which 
the template spectra were multiplied, which, for normal- 
ized spectra, is the fractional contribution of the secondary 
star to the total light in eclipse. These results are also sum- 
marised in Tabled] 

For AC Cnc, the spectral types with the lowest x^ val- 
ues are G7V and K2V in blue and G9V in red, by no means 
offering a definitive answer. However, the fractional contri- 
bution of the secondary star must be less than one, ruling 
out a G type companion. By visually inspecting each of the 
spectra, we settle on a spectral type for the secondary star 
in AC Cnc of K2 ± IV. A plot of the AC Cnc average eclipse 
spectrum, a broadened template spectrum and the residual 
of the optimal subtraction is shown in Fig. |7| The analy- 
sis using a K2V template results in a t; sin i measurement 
of 136 kms~^ in blue and 134 kms~^ in red, prompting us 
to adopt usinj — 135 ± 3 kms~^. This encompasses the 
vsini value for all the G and early-mid K templates (ex- 
cept for the blue G5V) within 2a. The error also reflects all 
of the other variations noted at the beginning of the pre- 
vious paragraph. [Schlegel, Kaitchuck fc Honeycutt (1984)] 
conclude that the secondary is a late G or early K star, 
not later than K3. We further limit this to an early 
K star, most likely K2V, agreeing with the studies of 
[Yamasaki, Okazaki fc Kitamura (1983)| and [Zhang (1987)| 
The results, however, conflict with the K5 estimate of 
[Shugarov (1981)] based on UBV colours. We find that, in 
eclipse, the secondary star in AC Cnc contributes 85 ± 5 



per cent of the total light in the blue and 74 ± 19 per cent 
in the red, assuming a K2 ± IV spectral type. 

For V363 Aur, the G7V template yields the lowest x^ 
value in blue, and the G5V proves the best in red. Unfor- 
tunately, we did not record spectra of the G6V and G7V 
templates in red, so we can only conclude from this analysis 
and by visual inspection that V363 Aur has a secondary of 
G7 ± 2V. The average vsini is 147 ± 5 kms~^ in the blue 
and 139 ± 5 kms~^ in the red. We therefore adopt a com- 
promise value of vsini — 143 ± 5 kms~^, encompassing all 
V sin i measurements for a G or early-mid K type secondary 
star. [Schlegel, Honeycutt fc Kaitchuck (1986) [ conclude that 
the spectral type is late G, and can be no later than K3, in 
agreement with this study. We do, however, rule out the es- 
timate of a GOV star by Szkody fc Crosa (1981)[ Using our 
adopted spectral type of G7 ± 2V, we find that in eclipse 
the secondary contributes 25 ± 3 per cent of the light in the 
blue and 45 ± 8 per cent in the red. 

At first glance, the fractional contributions of the sec- 
ondary stars in the two systems during eclipse appear to 
be inconsistent. The K2 ± IV secondary star in AC Cnc 
contributes a larger fraction in blue than red, whereas the 
(intrinsically bluer) G7 ± 2V secondary in V363 Aur con- 
tributes a smaller fraction. This can be explained by con- 
sidering the diflerent geometries of the two systems. In AC 
Cnc, almost all of the disc is obscured during eclipse, leaving 
only the redder outer disc uneclipsed. In V363 Aur, however, 
a large portion of the blue inner parts of the disc are still 
visible during eclipse (seen in Fig. IllH . contributing signifi- 
cantly to the blue eclipse light. Outside eclipse, we measure 
the fractional contribution of the secondary star in AC Cnc 
to be 19 ± 2 per cent in blue and 40 ± 11 per cent in red. 
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Table 2. (a) Times of mid-eclipse for AC Cnc according to 
Kurochkin & Sliugarov (1980; KS80), Yamasaki, Okazaski & Ki- 
tamura (1983; YOK83), Schlegel, Kaitchuck & Honeycutt (1984; 
SKH84), Zhang (1987; Z87) and this paper, (b) Times of mid- 
eclipse for V363 Aur according to Home, Lanning &; Gomer (1982; 
HLG82), Schlegel, Honeycutt & Kaitchuck (1986; SHK86), Rut- 
ten, van Paradijs & Tinbergen (1992; RvPT92) and this paper. 
The uncertainties on the mid-eclipse times measured from our 
data are 0.001; all other uncertainties are taken to be 0.005. 

(a) AC Cnc 



Cycle 
(E) 



HJD 

at mid-eclipse 

(2,400,000+) 



O-C 

(sees) 



Reference 



-59456 


34059.348 


-37.03 


KS80 


-50866 


36640.446 


-336.27 


KS80 


-50813 


36656.367 


-708.29 


KS80 


-47342 


37699.327 


-474.74 


KS80 


-33860 


41750.380 


886.55 


KS80 


-31241 


42537.320 


-20.02 


KS80 


-30136 


42869.344 


-331.40 


KS80 


-29177 


43157.504 


-149.38 


KS80 


-29074 


43188.450 


-424.07 


KS80 


-29061 


43192.362 


76.43 


KS80 


-29038 


43199.275 


250.81 


KS80 


-28848 


43256.368 


447.87 


KS80 


-28785 


43275.293 


8.90 


KS80 


-28165 


43461.579 


-857.79 


KS80 


-26867 


43851.619 


891.34 


KS80 


-26598 


43932.439 


162.19 


KS80 


-26588 


43935.443 


95.25 


KS80 


-26578 


43938.448 


114.72 


KS80 


-26525 


43954.376 


347.49 


KS80 


-23098 


44984.1119 


313.96 


YOK83 


-23094 


44985.3141 


339.02 


YOK83 


-23088 


44987.1165 


298.86 


YOK83 


-21522 


45457.6637 


254.17 


SKH84 


-19342 


46112.7032 


134.62 


Z87 


-19335 


46114.8053 


27.29 


Z87 


-19312 


46121.7177 


149.82 


Z87 


-6 


51922.7339 


-8.03 


This Paper 


-3 


51923.6350 


-36.75 


This Paper 





51924.5367 


-13.63 


This Paper 



(b) V363 Aur 






44557.9495 


-168.03 


HLG82 


3 


44558.9128 


-204.81 


HLG82 


6 


44559.8772 


-146.54 


HLG82 


105 


44591.6813 


-46.80 


HLG82 


106 


44592.0023 


-67.70 


HLG82 


4694 


46065.8614 


51.42 


SHK86 


4700 


46067.7877 


-48.05 


SHK86 


4703 


46068.7544 


208.94 


SHK86 


4706 


46069.7182 


215.37 


SHK86 


9365 


47566.3817 


8.49 


RvPT92 


9368 


47567.3466 


109.96 


RvPT92 


10057 


47788.6818 


70.93 


RvPT92 


10060 


47789.6457 


85.99 


RvPT92 


10082 


47796.7142 


187.84 


RvPT92 


10088 


47798.6401 


53.81 


RvPT92 


22916 


51919.5301 


-12.64 


This Paper 


22922 


51921.4577 


0.21 


This Paper 



Similar values are found for V363 Aur with out-of-eclipse 
contributions of 19 ± 4 per cent in blue and 35 ± 7 per cent 
in red. 



3.8 The Kr Correction 

The irradiation of the secondary stars in CVs by the 
emission regions around the white dwarf and the bright 
spot has been shown to influence the measured Kr (e.g. 
I Wade fc Home 19881 and IWatson fc Dhillon 200Tll . For ex- 
ample, if absorption lines are quenched on the irradiated 
side of the secondary, the centre of Ught will be shifted to- 
wards the back of the star. The measured Kr will then be 
larger than the true (dynamical) value. 

We must now determine whether the secondary stars 
in AC Cnc and V363 Aur are irradiated, which can be ob- 
servationally tested in two ways. Firstly, the rotationally 
broadened line profile would be distorted if there was a non- 
uniform absorption distribution across the surface of the sec- 
ondary star | |Davey fc Sm ith 1992). This would result in a 
non-sinusoidal radial velocity curve. Secondly, one would ex- 
pect a depletion of absorption line flux from the secondary 
star at phase 0.5, where the quenched inner-hemisphere is 
pointed towards the observer (e.g. IFriend e t al. 1990}. We 
applied these tests to the AC Cnc and V363 Aur data. 

The secondary star radial velocity curves were pro- 
duced by cross-correlating the CV spectra with the best- 
fitting smeared and broadened template spectra as described 
in Section 13.61 This time the cross-correlation peaks were 
plotted against phase to produce the radial velocity curves 
shown in the lower panels of Fig. |H] The radial velocity 
curves of both AC Cnc and V363 Aur are clearly eccen- 
tric in comparison to the sinusoidal fits represented by the 
thin solid lines. 

The variation of secondary star absorption line fiux with 
phase for AC Cnc and V363 Aur is shown in the top pan- 
els of Fig. |5] These lightcurves were produced by optimally 
subtracting the smeared and rotationally broadened best- 
fitting template from the individual CV spectra (with the 
secondary radial velocity shifted out) as described in Sec- 
tion 13.71 This time, however, the spectra were continuum 
subtracted rather than normalised to ensure that the mea- 
surements were not affected by a fluctuating disc brightness. 
The constants produced by the optimal subtraction are sec- 
ondary star absorption line fluxes, correct relative to each 
other, but not in an absolute sense. The dashed lines super- 
imposed on the lightcurves represent the variation of flux 
with phase for a Roche lobe with a uniform absorption dis- 
tribution. The sinusoidal nature is the result of the chang- 
ing projected area of the Roche lobe through the orbit. The 
lightcurves of AC Cnc and V363 Aur exhibit a drop in flux 
at phase 0.5 in comparison with the uniform Roche lobe. 

These two pieces of evidence, as well as the observed 
Balmer emission from the inner hemisphere of the secondary 
stars seen in Figs.j^and |1| and the weakening of the CCFs 
around phase 0.5 seen in Fig. |S| suggest that the secondary 
stars in AC Cnc and V363 Aur are indeed irradiated and we 
must correct the Kr values accordingly. 

It is possible to correct Kr for the effects of irradia- 
tion by modelling the secondary star absorption line flux 
distribution. In our model, we divided the secondary Roche 
lobe into 40 vertical slices of equal width. We then pro- 
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Figure 7. Orbitally-corrected average eclipse spectra of the CVs (top) with the best— fitting broadened template (middle) and the 
residuals after optimal subtraction (bottom). The template spectra have been multiplied by the optimal factor found from the optimal 
subtraction. All of the spectra are normalised and have been offset on the plots by an arbitrary amount for clarity. Left: AC Cnc 
(eclipse phases —0.05 < (f> < 0.05 in blue and red) with a K2V template. Right: V363 Aur (eclipse phases —0.07 < ij> < 0.09 in blue 
and —0.09 < <j> < 0.09 in red) with a G5V (red) and G7V (blue) template. The wavelength limits of the plots are those used for the 
cross— correlation and optimal extraction procedures for each object. The region between the dotted lines in the red spectrum of V363 Aur 
was omitted from the cross-correlation and optimal subtraction procedures due to the presence of a weak emission line. The average 
eclipse spectra shown (and the corresponding phase ranges listed above) are those used to determine the fractional contribution of the 
secondary stars given in Table l4l 



duced a series of model lightcurves, varying the numbers of 
slices omitted from the inner hemisphere of the secondary 
which contribute to the total flux. The model lightcurves 
were then scaled to match the observed data, and the 
best-fitting model found by measuring the y^ between the 
two. In all models, we used a gravity-darkening parame- 
ter j3 = 0.08 and limb-darkening coefficient u = 0.5 (e.g. 
IWatson fc Dhillon 200l|l . There is evidence for a secondary 
eclipse in both systems, so we have omitted points around 
phase 0.5 from the fits. (We tried a model which included 
an accretion disc to reproduce the secondary eclipse, but the 
results were exactly the same as omitting the points.) Once 
the best-fitting lightcurve was found, we produced fake CV 
spectra from the models, which were cross-correlated with 
a fake template star to produce a synthetic radial velocity 
curve. In the first instance, the synthetic curve mimicked 
the non-sinusoidal nature of the observed data, but with 
a larger semi-amplitude. This was expected, as the model 
input parameters used the uncorrected Kr derived in Sec- 
tion 13.101 We then lowered Kr and repeated the process, 
until the semi-amplitude of the model and observed radial 
velocity curves were in agreement, each time checking the 
lightcurve models for goodness of fit. The resulting Kr was 
then adopted as the real (or dynamical) Kr value. 

For AC Cnc, the best-fitting model lightcurve was pro- 
duced by omitting 8 slices when fitting both the blue and 
red data. The model lightcurves omitting 7, 8 and 9 slices 
are shown by the solid lines in Fig.|H| Our final model, which 
has an input Kr of 176 kms~^, produces the radial velocity 
curve shown as the thick solid line in Fig. |H| There is ex- 
cellent agreement between this and the observed data. The 
best-fitting lightcurve models for V363 Aur have 10 slices 
omitted in blue and 8 slices in red. Model lightcurves omit- 
ting 9, 10 and 11 slices in blue and 7, 8 and 9 slices in red 



are once again shown as solid lines in Fig. |H] The red data 
are very noisy, so we used the blue data to obtain a cor- 
rected Kr of 168 kms~^. The model with this input value 
has the radial velocity curve plotted as the thick solid line 
in Fig. |H1 It should be noted that if gravity-darkening and 
limb-darkening are neglected, the best fit lightcurves remain 
the same in all cases, and produce Kr values which are ~ 3 
kms~^ lower. 

The points around primary eclipse in the blue lightcurve 
of AC Cnc were also omitted from the above fit, as they show 
a very sharp decrease in flux. Although some of this can be 
attributed to the reduced projected area of the secondary 
at phase 0, the feature is too sharp and too deep for this 
to be the only explanation. It is likely that this feature is 
an artefact of the slit -loss correction procedure, where B- 
band photometry has been used to correct secondary fea- 
tures which actually lie closer to the V-band. Because the 
B-band eclipse is deeper than the V-band eclipse, we get 
a residual sharp dip in the secondary star flux at phase 0. 
There is no corresponding dip in the red data, as the pho- 
tometry and spectroscopy wavelengths are closely matched. 
There is no corresponding dip in flux at phase in the blue 
V363 Aur data, even though both objects have been re- 
duced in the same way. We believe this is because the V- 
band eclipse depth is approximately the same as the B-band 
eclipse depth in V363 Aur. 

In summary, we correct the Kr of AC Cnc from 186 
kms~^ to 176 kms"^ and the Kr of V363 Aur from 184 
kms~^ to 168 kms~^. 



3.9 The distances to AC Cnc and V363 Aur 

By finding the apparent magnitude of the secondary star 
from its contribution to the total light during eclipse, and 
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Figure 8. Upper panel: Secondary star lightcurves with model fits (solid lines). For AC Cnc, model fits are shown for Roche lobes with 
7, 8 and 9 slices removed in both blue (left) and red (right). For V363 Aur, model fits are shown for Roche lobes with 9, 10 and 11 slices 
removed in blue (left) and 7, 8 and 9 in red (right). The lower the line, the more slices removed. The dashed line in all plots represents 
a model where slices have been removed. The data have been phase— binned into 50 bins to increase signal-to-noise. Lower panel: 
Measured secondary star radial— velocity curve with a sinusoidal fit (thin solid line) and the best-fitting model fit (thick solid line). In 
all panels, the open circles indicate points that have been omitted from the fits and all data have been folded to show 2 orbital phases. 
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Table 4. vsini values (and miniinuin reduced x^) for AC Cnc and V363 Aur cross-correlated with the rotationally-broadened profiles 
of G5-K8 dwarf templates. Degrees of freedom: AC Cnc, 1143 (blue), 498 (red); V363 Aur, 580 (blue), 360 (red). Also shown is the 
fractional contribution of the secondary star to the total light during the eclipse phases, and the position of the strongest peak in the 
skewinaps when cross— correlated with each template using a 7-velocity of 40 kms~^ for AC Cnc and kms~^ for V363 Aur. 

(a) AC Cnc 



Templates 


min x^ 


V sin j 


min x^ 


f sin j 


Fractional 


Fractional 


(Kx,Ky) from 


{K^,Ky) from 




(blue) 


at min x 


(red) 


at min x^ 


contribution 


contribution 


skew map 


skew map 






(blue) 




(red) 


of secondary 


of secondary 


(blue) 


(red) 






kms-l 




kms-i 


(blue) 


(red) 


kms^-"- 


kms -•- 


G5V 


1.081 


142 


1.106 


136 


1.24 ±0.03 


1.13 ±0.10 


(3,187) 


(-4,185) 


G6V 


1.066 


140 


- 


- 


1.15 ±0.03 


- 


(3,186) 


- 


G7V 


1.038 


140 


- 


- 


1.18 ±0.03 


- 


(-2,185) 


- 


G8V 


1.046 


141 


1.100 


137 


1.20 ±0.03 


1.22 ±0.11 


(-1,186) 


(-4,187) 


G9V 


1.080 


139 


1.097 


136 


1.19 ±0.03 


1.37 ±0.12 


(3,186) 


(-4,186) 


KOV 


1.043 


138 


1.107 


135 


1.04 ±0.02 


0.99 ±0.09 


(1,185) 


(-4,187) 


KIV 


- 


- 


1.112 


134 


- 


0.93 ±0.08 


- 


(-9,186) 


K2V 


1.038 


136 


1.112 


134 


0.85 ± 0.02 


0.74 ±0.07 


(0,187) 


(0,185) 


K3V 


1.127 


134 


- 


- 


0.82 ± 0.02 


- 


(-2,186) 


- 


K4V 


1.406 


138 


1.137 


131 


0.62 ± 0.02 


0.46 ±0.04 


(3,189) 


(-2,185) 


K5V 


1.288 


135 


1.128 


131 


0.67 ±0.02 


0.56 ±0.05 


(0,187) 


(-10,186) 


K7V 


1.319 


136 


1.125 


130 


0.65 ± 0.02 


0.55 ±0.05 


(4,186) 


(-3,186) 


K8V 


1.437 


136 


1.138 


129 


0.63 ± 0.02 


0.50 ±0.05 


(0,188) 


(-4,184) 


(b) V363 Aur 


G5V 


1.081 


148 


1.119 


139 


0.25 ±0.01 


0.43 ±0.07 


(-7,184) 


(10,185) 


G6V 


1.081 


146 


- 


- 


0.25 ±0.01 


- 


(-8,184) 


- 


G7V 


1.054 


148 


- 


- 


0.27 ±0.01 


- 


(-13,181) 


- 


G8V 


1.068 


148 


1.122 


138 


0.26 ±0.01 


0.46 ±0.08 


(-11,184) 


(12,184) 


G9V 


1.098 


143 


1.120 


140 


0.24 ±0.01 


0.52 ±0.09 


(-4,183) 


(13,185) 


KOV 


1.086 


143 


1.126 


140 


0.21 ±0.01 


0.37 ±0.07 


(-6,185) 


(15,184) 


KIV 


- 


- 


1.120 


138 


- 


0.37 ±0.06 


- 


(7,189) 


K2V 


1.107 


142 


1.136 


146 


0.17 ±0.01 


0.29 ±0.06 


(-9,186) 


(13,185) 


K3V 


1.150 


139 


- 


- 


0.17 ±0.01 


- 


(-10,184) 


- 


K4V 


1.294 


134 


1.140 


146 


0.10 ±0.01 


0.19 ±0.04 


(0,191) 


(9,186) 


K5V 


1.234 


134 


1.132 


141 


0.11 ±0.01 


0.24 ±0.04 


(5,188) 


(8,188) 


K7V 


1.252 


132 


1.134 


145 


0.11 ±0.01 


0.23 ±0.04 


(0,191) 


(12,187) 


K8V 


1.297 


134 


1.138 


148 


0.10 ±0.01 


0.22 ±0.04 


(-4,188) 


(8,187) 



estimating its absolute magnitude, we can calculate the dis- 
tance (d) to each system using the equation: 



51og(d/10) =mv -Mv -d Av/WOO 



(5) 



where Ay is the visual interstellar extinction in magnitudes 
per kpc. 

There are a number of ways of estimating the absolute 
magnitude of the secondary star, assuming it is on the main 
sequence (e.g. IPatterson 19841 IWarner 1995bl [Gray 1992| l. 
The distance estimates given below take into account all of 
these techniques. 

Another method of finding the distance is to determine 
the angular diameter of the secondary star from the observed 
flux and a surface brightness calibration that we derive from 
the Barnes-Evans relation (Barnes & Evans 197^, 

F„ = 4.2207 - 0.1 M) - 0.5 log (f) = 3.977 - 0.429( V - -R)o (6) 

where Vo and {V — R)o are the unreddened V magnitude 
and {V — R) colour of the secondary star, and </> is the stellar 
angular diameter in arc milliseconds. 



3.9.1 AC Cnc 

During the eclipse phases given in Section [3.71 the average 
apparent magnitude of AC Cnc is R — 14.7 ± 0.1, of which 
the secondary contributes 74 ± 19 per cent, and B = 15.9 ± 
0.1, of which the secondary star contributes 85 ± 5 per cent. 
The apparent magnitude of the secondary is therefore R — 
15.0 ± 0.3 and B = 16.1 ± 0.1. Using typical ■B-F and V~R 
values for an early K star from |Gray (1992)] we arrive at an 
apparent V magnitude of 15.5 ± 0.3. We adopt an absolute 
magnitude of My = +6.8 ± 0.5 as an average of the various 
methods referenced in the previous paragraph. Assuming 
zero interstellar extinction Ilia Dous 199111 . the distance to 
AC Cnc is 550 ± 150 pc. 

Using the Barnes-Evans relation with a (V — _R)o value 
typical of an early K star (0.74 ± 0.10; [Gray 1992^ , the 
Vo value of 15.5 ± 0.3 found above and the radius of the 
secondary star derived in Section r3.10l we obtain a distance 
of 750 ± 250 pc. 

Published estimates of the distance to AC Cnc 
agree with our calculated values. [Shugarov (1981)] calcu- 
lates a distance of 480 pc assuming a K5 secondary 
star. [Patterson (1984)] derives a value of 400 pc by com- 
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billing an Mv-H/3 equivalent-width relationship, prop- 
erties of the secondary and the continuum shape of 
the spectrum. [Warner (1987)| calculates a distance of 
800 pc using the secondary star characteristics from 
[Yamasaki, Okazaki fc Kitamura (1983)| and [Zhang (1987)] 
concludes d = 500 ± 100 pc, again estimating the secondary 
star properties. 

3.9.2 V363 Aur 

The average apparent magnitude of V363 Aur during the 
eclipse phases quoted in Section [3.71 is B = 14.8 ±0.1 and 
R = 13.9 ± 0.1, of which the G7 ± 2V secondary star 
contributes 25 ± 3 per cent and 45 ± 8 per cent, respec- 
tively. This corresponds to an apparent magnitude of the 
secondary star of B = 16.3 ± 0.2, R = 14.8 ± 0.2 and, 
assuming typical B — V and V — R values for a G7 ± 
2V star | |Gray 1992^ , we calculate V = 15.4 ± 0.2. There 
is evidence of interstellar absorption in the average spec- 
trum (Fig. 0, suggesting that it is important to take ex- 
tinction into account. [Szkody fc Crosa (1981)] give a value 
of E{B — V) = 0.3, although the UV spectrum used was un- 
derexposed and noisy. [la Dous (1991)[ however, finds no ex- 
tinction to V363 Aur, a result contradicted by measurements 
of [Rutten, van Paradijs fc Tinbergen (1992)[ who measure 
E{B — V) = 0.1. We adopt this value, which results in 
Av = 0.32 IjSchefHer 1982|l . We use Mv = +6.0 ± 0.5 as 
an average absolute magnitude of the G7 ± 2V star sec- 
ondary star in V363 Aur, which results in a distance of 700 
± 250 pc. 

Using the Barnes-Evans relation with a {V — R)o value 
typical of a mid-late G star (0.56 ± 0.1; [Gray 1992^ , a Vo 
value of 15.1 ± 0.2 assuming extinction to be 0.3 mag and 
the radius of the secondary star derived in Section f3.10l we 
obtain a distance of 1000 ± 250 pc. 

[Szkody fc Crosa (1981)[ estimate a distance to 
V363 Aur of 900 pc assuming a GO dwarf secondary. 
[Patterson (1984)] derives an uncertain value of 1000 pc 
from the My-H/3 equivalent -width relationship, interstellar 
absorption and the continuum shape of the spectrum. In his 
compilation of distances, [Berriman (1987)] gives a distance 
of 1100 pc from disc properties, 900-1300 pc from red star 
spectrophotometry and >450pc from the infrared properties 
of the secondary. [Rutten, van Paradijs fc Tinbergen (1992)] 
calculate a distance of 530 pc from a black body fit to the 
spectrum of the central part of the disc and 720 pc using a 
value for the fractional contribution of the secondary. 

3.10 System Parameters 

Using the Kr and v sin i values found in Sections l3.7l and l3.8l 
in conjunction with the period determined in Section 13.11 
and a measurement of the eclipse full-width at half depth 
{A(j)i/2), we can calculate accurate system parameters for 
AC Cnc and V363 Aur. 

In order to determine A(;/>i/2, we estimated the flux 
out of eclipse (the principal source of error) and at eclipse 
minimum, and then measured the full-width of the eclipse 
half-way between these points. The eclipse full-width at 
half-depth was measured to be A</)i/2 = 0.096 ± 0.002 for 
AC Cnc and 0.063 ± 0.002 for V363 Aur from the B and 
R-hand lightcurves in Fig. [5] 



We have opted for a Monte Carlo approach similar to 
[Home, Welsh fc Wade (1993)[ to calculate the system pa- 
rameters and their errors. For a given set of Kr, vsini, 
A(f)i/2 and P, the other system parameters are calculated 
as follows. 

R2/a can be estimated because we know that the sec- 
ondary star fills its Roche lobe (as there is an accretion disc 
present and hence mass transfer) . R2 is the equatorial radius 
of the secondary star and a is the binary separation. We used 
Eggleton's formula ( [Eggleton 1983^ which gives the volume- 
equivalent radius of the Roche lobe to better than 1 per cent, 
which is close to the equatorial radius of the secondary star 
as seen during eclipse. 



0.49g^''^ 



(7) 



a 0.6g2/3-f ln(l + gi/3 

The secondary star rotates synchronously with the orbital 
motion, so we can combine Kr and ?;sini, to get 

^(1 + .) = ^. (8) 

a Kr 

By considering the geometry of a point eclipse by a spherical 
body (e.g. [Dhillon, Marsh fc Jones 1991| l, the radius of the 
secondary can be shown to be 



— I — sin 7rA(^i/2 -I- cos nAcj>i/2 cos i, 
a J 



(9) 



which, using the value of R2/a obtained using equations [7] 
and [HI allows us to calculate the inclination (i) of the sys- 
tem. The geometry of a disc eclipse can be approximated 
to a point eclipse if the light distribution around the white 
dwarf is axi-symmetric (e.g. lDhillon 1990ll . This approxima- 
tion is justified given the symmetry of the primary eclipses 
in the photometry light curves (Figure [2J . Kepler's Third 
Law gives us 



K\Porb _ Misin^i 
~^i^G~' (1 + g)" 



(10) 



which, with the values of q and i calculated using equations 
E|[H]and[5] gives the mass of the primary star. The mass of 
the secondary star can then be obtained using 



M2 = qMi 



(11) 



The radius of the secondary star is obtained from the equa- 
tion 

II sin i 2-K sin i 



R2 P ' (^2) 

(e.g. IWarner 1995afl and the separation of the components, 
a, is calculated from equations |S]and I12l with q and i now 
known. 

The Monte Carlo simulation takes 10 000 values of Kr, 
II sin i, and A(^i/2 (the error on the period is deemed to be 
negligible in comparison to the errors on Kr, usini, and 
A4'i/2), treating each as being normally distributed about 
their measured values with standard deviations equal to the 
errors on the measurements. We then calculate the masses of 
the components, the inclination of the system, the radius of 
the secondary star, and the separation of the components, as 
outlined above, omitting {Kr, vsini, A(/)i/2) triplets which 
are inconsistent with sin i < 1 . Each accepted Mi , M2 pair is 
then plotted as a point in Figure|^ and the masses and their 
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errors are computed from the mean and standard deviation 
of the distribution of these pairs. 

In the case of AC Cnc, we find that Mi = O.76±O.O3M0 
and Kh = 0.77 ± O.OSMq; for V363 Aur we find Afi = 
0.90 ± Q.O6M0 and Ma = 1.06 ± O.IIMq. The values of 
all the system parameters deduced from the Monte Carlo 
computation are listed in Table |K| including ifij-corrected 
and non i^jj-corrected values for comparison. 

We computed the radius of the accretion discs in 
AC Cnc and V363 Aur using the geometric method outlined 
in |Dhillon, Marsh fc Jones (1991)| The phase half- width of 
eclipse at maximum intensity was found in Section r3.3l to be 
A0 = 0.09 ± 0.01 for AC Cnc and Ac/) = 0.078 ± 0.005 for 
V363 Aur. Combining A(j) with q and i derived above pro- 
duces an accretion disc radius (Rd) in terms of the volume 
radius of the primary's Roche lobe (-Ri). We find accretion 
disc radii of Rd/Ri = 0.61 ± 0.14 and Rd/Ri = 0.71 ± 0.17 
for AC Cnc and V363 Aur, respectively, which are lower than 
those quoted by |Harrop-Allin fc Warner (1996)] but consis- 
tent within the errors; AC Cnc: 0.70 < Rd/Ri < 0.86, 
V363 Aur: Rd/Ri > 0.92. 

The empirical relation obtained by 

[Smith fc Dhillon (1998)| between mass and radius for 
the secondary stars in CVs is given by. 



A ^(0.93 ±0.09)^ 
Rq ^ \Mq ^ 



(0.06 ±0.03). 



(13) 



This predicts that if the secondary stars in AC Cnc and 
V363 Aur are on the main-sequence, they should have radii 
of O.78i?0 and l.O57?0, respectively. These values agree with 
our measured values of 0.83 ± 0.03 -R0 and 0.97 ± 0.04 Rq to 
within the errors. IGra y (1992)| gives M = O.76M0 and R = 
O.75-R0 for a K2 dwarf and M = O.98M0 and R = O.96i?0 
for a G5 dwarf, also in agreement with our measured values. 
We conclude that the secondaries in AC Cnc and V363 Aur 
are similar to main-sequence stars. 



4 DISCUSSION 

4.1 Are AC Cnc and V363 Aur SW Sex stars? 

The SW Sex stars are a sub-class of NLs which have peculiar 
spectral properties - see lHellier 20001 and references therein 
for a summary of the current models. Most SW Sex stars 
share the following properties (e.g. IHoard et al. 2003t : 

(i) They are usually high-inclination, eclipsing systems. 

(ii) Their spectra exhibit single-peaked emission lines 
rather than the double-peaked lines expected from near 
edge-on discs. 

(iii) The Balmer and Hel emission lines usually contain 
transient absorption cores, especially around photometric 
phase 0.5. 

(iv) They have high levels of excitation, with Hell A4686A 
emission often comparable in strength to H/3. 

(v) The low-excitation lines (Balmer and Hel) exhibit 
shallow or absent eclipses. 

(vi) High velocity S-waves are often seen in the trailed 
spectra. 



(vii) The emission-line radial velocity curves show large 
phase shifts between spectroscopic conjunction and photo- 
metric mid-eclipse. 

AC Cnc and V363 Aur show many of the features de- 
scribed above. V363 Aur shows all of the features, so must be 
considered a definite SW Sex star. AC Cnc does show high 
excitation features, but not to the extent of V363 Aur. Tran- 
sient absorption features are seen in low-excitation lines of 
both systems. An interesting absorption feature also occurs 
around phase in V363 Aur, which deserves further discus- 
sion. 

The low-excitation emission lines of V363 Aur show the 
broad primary eclipse expected from a gradual obscuration 
of the accretion disc by the secondary (see Figl^J. Around 
phase 0, however, there is a sharp increase in Balmer line 
emission followed by a rapid decrease. The feature comes and 
goes too rapidly to be attributable to a permanent emission 
source on the back of the secondary star. We suggest that 
this increase in flux can be explained by the eclipse of an op- 
tically thick region at the centre of the disc. The eclipse of 
this region, which would have an absorption line spectrum, 
would result in an overall increase in the line flux at the 
phases observed. Fig. llUl is the trailed spectrum of H/3 mag- 
nified to show the effect more clearly. The blue wing of the 
line increases in flux first, in agreement with an eclipse of the 
centre of the disc. Fig. 1111 is a reconstruction of V363 Aur 
using the isTij-corrected system parameters listed in Fig. |21 
The left-hand diagram shows the system just before the in- 
crease in flux at phase 0.97, where the centre of the disc is 
just about to be eclipsed by the secondary. In the central dia- 
gram, the centre of the disc is eclipsed, corresponding to the 
brief rise in line flux. The right-hand diagram shows the sys- 
tem at phase 1.03, when the centre of the disc comes out of 
eclipse, and the line flux drops again. The short timescale of 
this feature is explained by the grazing nature of the eclipse 
of the optically thick region at the centre of the disc. This 
model could also account for the behaviour of other SW Sex 
stars where the Balmer lines show only a shallow eclipse 
com pared to the continuum. [Dhillon, Marsh fc Jones (1997)| 
and [Groot, Rutten fc van Paradijs (2001)[ suggest a similar 
explanation for SW Sex. The latter authors invoked an opti- 
cally thick absorption-line source coincident with the bright 
spot, similar to that of a late B or early A-type star. 



4.2 Mass Transfer Stability 

The mass ratio (q = M2/M1) of a CV is of great signifi- 
cance, as it governs the properties of mass transfer from the 
secondary to the white dwarf primary. This in turn governs 
the evolution and behaviour of the system. 

The secondary star has two timescales on which it re- 
sponds to mass loss. Firstly, the star returns to hydrostatic 
equilibrium on a dynamical timescale, which is the sound- 
crossing time of the region affected. Secondly, on a longer 
timescale, the star settles into a new thermal equilibrium 
configuration. The second timescale a star responds on is 
therefore the thermal, or Kelvin-Helmholtz, timescale. 

The two timescales upon which the secondary responds 
to mass loss leads to two types of mass transfer instability. 
If, upon mass loss, the dynamical response of the secondary 
is to expand relative to the Roche lobe, mass transfer is 
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Figure 9. Monte Carlo determination of system paramters for AC Cnc (left) and V363 Aur (right). Each dot represents an Mi,M2 
pair; the solid curves satisfy the i)sini and corrected Kn constraints, and the dashed lines mark lines of constant inclination. 

Table 5. System parameters for AC Cnc and V363 Aur. The Monte Carlo results for corrected and uncorrected Kr values are shown 
for comparison. The radial velocity of the white dwarf (Kw) has also been calculated from the secondary star parameters. Distances 
are quoted using both techniques described in the text; by estimating the secondary star magnitude {d"'^'^) and using the Barnes- Evans 
relation (d^-^). 



Parameter 




AC Cnc 






V363 Aur 






Measured 


Monte Carlo 


Kfl-corrected 


Measured 


Monte Carlo 


iiT/j -corrected 




Value 


Value 


Value 


Value 


Value 


Value 


Porb (d) 


0.30047747 




0.30047747 


0.32124187 




0.32124187 


Kr (kms-1) 


186 ±3 


186 ±3 


176 ±3 


184 ±5 


184 ±5 


168 ±5 


vsmi (kms"-*^) 


135 ±3 


135 ±3 


135 ±3 


143 ±5 


143 ±5 


143 ±5 


^^01/2 


0.096 ± 0.003 


0.096 ± 0.003 


0.096 ± 0.003 


0.063 ± 0.002 


0.063 ± 0.002 


0.063 ± 0.002 


1 




0.94 ±0.04 


1.02 ±0.04 




1.04 ± 0.06 


1.17 ±0.07 


i° 




76.3 ±0.8 


75.6 ±0.7 




70.5 ± 0.4 


69.7 ±0.4 


Kyv (kms-1) 




175 ±6 


179 ±6 




190 ±9 


196 ±9 


Mi/Mq 




0.82 ± 0.04 


0.76 ±0.03 




1.03 ±0.07 


0.90 ±0.06 


M2/Mq 




0.78 ± 0.05 


0.77 ±0.05 




1.06 ±0.11 


1.06 ±0.11 


R2/Rq 




0.82 ±0.02 


0.83 ±0.02 




0.96 ±0.04 


0.97 ±0.04 


a/Re 




2.21 ±0.04 


2.18 ±0.04 




2.52 ±0.07 


2.47 ±0.07 


d""" (pc) 


550 ±150 






700 ± 250 






d^-E (pc) 


750 ± 250 






1000 ± 250 






Spectral type 


K2±l V 






G7±2 V 






of secondary 














Afli 


0.09 ±0.01 






0.078 ± 0.005 






Rd/Ri 






0.61 ±0.14 






0.71 ±0.17 



dynamically unstable and mass transfer proceeds on the dy- 
namical times cale. [Politano (1996)] made an analytic fit to 
the models of |Hjellming (1989)| to give the limit of dynam- 
ically stable mass loss, plotted as the solid line in Fig. 1121 
Dynamically stable mass transfer can occur if the CV lies be- 
low this line. This limit is important for low mass secondary 
stars (M2 < O.SA/q), as they have significant convective 
envelopes that tend to expand adiabatically in response to 
mass loss llde Kool 19921 . 

Thermally unstable mass transfer is possible if the dy- 
namic response of the star to mass loss is to shrink rela- 
tive to its Roche lobe (i.e. mass transfer is dynamically sta- 



ble). This occurs at high donor masses (A/2 > O.8M0) when 
the star has a negligible convective envelope and its adia- 
batic response to mass loss is to shrink (e.g. Ide Kool 19921 
IPolitano 199611 . Mass transfer then initially breaks contact 
and the star begins to settle into its new thermal equilibrium 
configuration. If the star's thermal equilibrium radius is now 
bigger than the Roche lobe, mass transfer is again unstable, 
but proceeds on the slower, thermal timescale. The limit of 
thermally stable mass transfer can be found by differenti- 
ating the main-sequence mass-radius relationship given in 
|Politano(1996), Thermally stable mass transfer can occur 
if the CV appears below the dotted line plotted in Fig. 1121 
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Figure 11. A model of V363 Aur using the i^ij-corrected system parameters listed in Table l5l During primary eclipse, the secondary 
star grazes the absorption region at the centre of the disc, resulting in the rapid increase and then rapid decrease in Balmer flux seen 
in the trailed spectra at phase 0. The left-hand diagram shows the system at phase 0.97, the centre diagram at phase 1.00 and the 
right-hand diagram at phase 1.03. 




Smith & Dhillon 199B 

Smith & Dhillon 199B [Kj only) 
ir North et al. 3000 
O Watson et al. 2003 

This paper 



Velocity (km/s) 

Figure 10. The V363 Aur H/3 trailed spectrum, magnified to 
show the phases 0.8-1.2. The emission feature is clearly seen 
around phase 0, moving from blue to red. 



Most CVs on the plot fall below both curves, implying 
that mass transfer is dynamically and thermally stable, as 
expected (the mass transfer rates observed in CVs are too 
low for unstable mass transfer to be occurring). An excep- 
tion is DX And, which appears to be dynamically unstable. 
This is, of course, not possible; a system undergoing dynam- 
ical mass transfer would rapidly form a common envelope. 
A solution is found in the fact that DX And has an evolved 
secondary star ( |Drew, Jones fc Woods 1993| ; the dynami- 
cal and thermal solutions plotted on b'\e. Vl'A are for ZAMS 
stars, so cannot be applied to evolved stars. 

The previous study of the masses of AC Cnc by 
[Schlegel, Kaitchuck fc Honeycutt (1984)] found a q value of 
1.24 ± 0.08, close to the limit for thermally stable mass 
transfer. In this study, we find a lower value of g = 
1.02 ± 0.04, placing it well within the theoretical limit. 
In the case of V363 Aur, we find a higher q value of 
1.1 7 ± 0.07 t han the value of g = 0.90 ± 0.10 quoted 
by [SchlegelTH oncycut t fc Kaitchuck (1986)] but it is still 
within the stability limit. The mass ratios found in this pa- 
per therefore place AC Cnc and V363 Aur within the region 
allowed by the theoretical constraints for stable mass trans- 
fer. 




-0.5 0.5 

log 10 (Mg/Mg) 

Figure 12. Critical mass ratios for mass transfer stability. The 
dotted line represents the condition for thermal instability; the 
solid line represents the condition for dynamical instability (Poli- 
tano 1996). Both curves assume the star is initially in ther- 
mal equilibrium. Mass ratios and secondary masses from the 
compilation of [Smith fc Dh illon (1998) North ct al . (2000) | and 
[Watson et al. (2003)| are overplotted. The mass ratios and sec- 
ondary star masses of AC Cnc and V363 Aur determined in this 
paper are also plotted. 
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